Multi-quark hadrons from Heavy Ion Collisions 
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Identifying hadronic molecular states and/or hadrons with multi-quark components either with 
or without exotic quantum numbers is a long standing challenge in hadronic physics. We suggest 
that studying the production of these hadrons in relativistic heavy ion collisions offer a promising 
resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their 
structures. Using the coalescence model for hadron production, we find that compared to the case 
of a non-exotic hadron with normal quark numbers, the yield of an exotic hadron is typically an 
order of magnitude smaller when it is a compact multi-quark state and a factor of two or more 
larger when it is a loosely bound hadronic molecule. We further find that due to the appreciable 
numbers of charm and bottom quarks produced in heavy ion collisions at RHIC and even larger 
numbers expected at LHC, some of the newly proposed heavy exotic states could be produced and 
realistically measured in these experiments. 

PACS numbers: 14.40.Rt,24.10.Pa,25.75.Dw 



(N 
> 
(N 

in 
oo 
o 



Finding hadrons with configurations other than the 
usual qq configuration for a meson and qqq for a baryon is 
a long standing challenge in hadronic physics. In 1970's, 
the tetraquark picture jl| was suggested as an attempt 
to understand the inverted mass spectrum of the scalar 
nonet. At the same time, the exotic H dibaryon Q was 
proposed on the basis of the color-spin interaction. While 
results from the long search for the H dibaryon in various 
experiments turned out to be negative, we are witnessing 
a renewed interest in this subject as the properties of sev- 
eral newly observed heavy states, including D s j(2317) Q 
and X(3872) cannot be properly explained within the 
simple quark model. 

An important aspect in understanding a multi-quark 
hadron involves the discrimination between a compact 
multi-quark configuration and a loosely bound molecular 
configuration with or without exotic quantum numbers. 
While the wave function of a loosely bound molecular 
configuration is dominantly composed of a bound state 
of well separated hadrons, the main Fock component of 
a compact multi-quark configuration typically has the 
size of a hadron, with little if any separable color sin- 
glet components. For a crypto-exotic state, one further 
has to distinguish it from a normal quark configuration. 
For example, /o(980) and ao(980) could be either normal 
quark-antiquark states [H , compact tetra-quark states [l| 
or weakly bound KK molecules |f|. 



Previously, discriminating between different configu- 
rations for a hadron relied on information about the de- 
tailedproperties of the hadron and its decay or reaction 
rate [7J. Moreover, searches for exotic hadrons have usu- 
ally been pursued in reactions between elementary par- 
ticles. In this letter, we will show that measurements 
from heavy ion collisions at ultrarelativistic energies can 
provide new insights into the problem and give answers 
to some of the fundamental questions raised above 0- 
Toj l . In particular, we focus on the yields of multi-quark 
hadrons in heavy ion collisions. To carry out this task, 
we first use the statistical model [lH , which assumes that 
the produced matter in relativistic heavy ion collisions 
is in thcrmodynamical equilibrium and is known to de- 
scribe the relative yields of normal hadrons very well, 
to normalize the expected yields. We then use the coa- 
lescence model [12 [. which is based on the sudden ap- 
proximation by calculating the overlap of the density 
matrix of the constituents in an emission source with 
the Wigner function of the produced particle, to take 
into account the effects of the inner structure of hadrons, 
such as angular momentum [13j and the multiplicity of 
quarks The coalescence model has been extensively 
used to study both light nucleus production in nuclear re- 
actions fl4j and hadron production from the quark-gluon 
plasma produced in relativistic heavy ion collisions Hi- 
18(| . In particular, it has successfully explained the ob- 
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served enhancement of baryon production in the inter- 
mediate transverse momentum region fl9L |20| and the 
quark number scaling of the elliptic flow of identified 
hadrons [2l|, [22| as well as the yield of recently discovered 
antihypertritons in heavy ion collisions at RHIC (23j . 

In the statistical model, the number of produced 
hadrons of a given type h is given by 
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with gh being the degeneracy of the hadron, and 
Vh (Tr) the volume (temperature) of the source when 
statistical production occurs. The fugacity is jh = 
jnc+n Ee (jiBB+v, B s)/T H ^ w h ere q anc j s are the baryon and 

strangeness numbers of the hadron with corresponding 
chemical potentials [Lb and /is, and n c (n E ) the number of 
(anti-)charm quarks. For central Au+Au (Pb+Pb) col- 
lisions at = 200 GeV (5.5 TcV) at RHIC (LHC), 
values for these parameters have been determined in 
Refs. 0, [2^ for particles in one unit of central rapidity 
in an expanding fire-cylinder model: Vh = 1908 (5152) 
fm 3 , T H = 175 MeV, fi s = 10 (0) MeV, and fi B = 20 (0) 
MeV. We fix j c = 6.40 (15.8) by requiring the expected 
total charm quark number N c = 3 (20) extracted from 
initial hard scattering at RHIC (LHC) to be equal to the 
sum of the yields of D, D* , D s and A c estimated in the 
statistical model. We note that all statistically produced 
hadrons from this fire-cylinder are essentially in the cen- 
tral unit rapidity. 

In the coalescence model, the number of hadrons of 
type h produced from the coalescence of n constituents, 
based on harmonic oscillator wave functions for the 
hadron internal structure, is given by 
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if we use the non-relativistic approximation, neglect the 
transverse flow of produced matter, and consider only 
the central unit rapidity as in Refs. d, Q. In Eq.([2]), gi 
is the degeneracy of the ith constituent and Ni its num- 
ber taken to be iV u =245 (662) and V s =150 (405), and 
V =1000 (2700)fm 3 for RHIC (LHC) @j; k is (1) for 
a s(p)-wave constituent; and = 1/ with u be- 
ing the oscillator frequency and /ii the reduced mass de- 
fined by = m^ +1 + (X)j=i m i)~ 1 with TO Ui d(m s )=300 
(500)MeV. Eq.([2]) is applicable if the number of con- 
stituent particles Nj is greater than one, which is the 
case even for charmed and bottom particles in central 
heavy ion collisions at RHIC and LHC. Furthermore, 
Eq. ([2]) shows that the addition of a s-wave or p-wave 
u/d-quark with i = 1 leads to the coalescence factor of 
about 0.360 or 0.093, respectively. Therefore, hadrons 
with more constituents are generally suppressed, and the 



p-wave coalescence is hindered with respect to the s-wave 
coalescence [l3| . 

In applying the coalescence model to multi-quark 
hadron production, we fix the oscillator frequencies by 
requiring the coalescence model to reproduce the ref- 
erence normal hadron yields in the statistical model. 
This leads to oj = 550 MeV for hadrons composed 
of light quarks. For hadrons composed of light and 
strange (charm) quarks, we fix oj s (uj c ) to reproduce the 
yields of A(1115) (A c (2286)) in the statistical model. For 
the A c (2286) yield, we include the feed-down contribu- 
tion according to N%$$ = #£$286) + A%f 2455) + 

7V s c a (2520) + °- 67 x ^1^(2625) • Fit ting this yield to that 
calculated in the coalescence model, we obtain uj c = 385 
MeV for m c = 1500 MeV. Similarly, we get lo s = 519 
MeV from the A(1115) yield after including the feed- 
down from the octet and decuplet states. 

The yields for weakly bound hadronic molecules are 
estimated using the coalescence of hadrons at the kinetic 
freezeout point (T F = 125 MeV, V F = 11322 (30569) fm 3 
for RHIC (LHC)). If the radius for hadronic molecules 
is known, the oscillator frequency uj can be fixed by 
uj = 3/(2/ii(r 2 )) for the 2-body s-wave state. If only 
the binding energy is given, we use the relation B.E. ~ 
Ti 2 / {2hiclq) and (r 2 ) ~ io/2: with ao being the s— wave 
scattering length, between the binding energy and the 
rms radius to obtain u> = 6 x B.E.. For example, 
for / (980), w /o(980) = 67.8 MeV using B.E. /o(980) = 

11.3 MeV. Table Usum- 
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marizes the parameters and possible decay modes for a 
selection of multi-quark candidates as well as proposed 
states KKN^, KNN^, DN, and DNN 

The yields of states listed in Table Q] are summarized 
in Table HT1 For example, possible configurations of the 
/o(980) could be an ss or a uu and dd state in addition to 
crypto-exotic configurations discussed before. For most 
of the states considered here, the coalescence yield from 
the compact multi-quark state is an order of magnitude 
smaller than that from the usual quark configuration 
as the coalescence of additional quarks are suppressed. 
Also, for the same hadronic state, the coalescence yield 
from the molecular configuration is similar to or larger 
than that from the statistical model prediction. The sim- 
ilarity in the yields from the statistical model and the co- 
alescence model prediction for a molecular configuration, 
despite the difference in the production temperatures Tc 
and Tp, can be attributed to the larger size of the molec- 
ular configuration forming at lower temperature but at a 
larger volume; hence the ratio of volumes uf /V is similar. 
The predicted appreciable yields of hadronic molecules in 
relativistic heavy ion collisions are in sharp contrast to 
those in high energy pp collisions, where molecular con- 
figurations with small binding energy are hard to pro- 
duce, particularly at high px [2^|. Our results do not 
change much if different forms of hadron wave functions 
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TABLE I: List of multi-quark states. For hadron molecules, the oscillator frequency WmoI. is fixed using the binding energy 
((B)) or the inter-hadron distance ((R)). The WmoI. for last two states is taken from corresponding two-body system ((T)). 



Pariticle 


m 

(MeV) 9 I J-k 


2q/3q/6q 


Aq/5q/8q 


Mol. 


k>Mol. 

(MeV) 


decay 
mode 


/o(aoU) 
a (980) 
-D s (2317) 
X(3872) 


98U 1 U U+ 

980 3 1 0+ 
2317 1 0+ 
3872 3 1+ 


qq [L — 1) 
qq (L = 1) 
cs (L = 1) 


qqss 
qqss 
qqcs 
qqcc 


K A 

KK 
DK 
DD* 


Q (Ti\ 
( .O^-D J 

67.8(B) 
273(B) 
3.6(B) 


7T7T (strong decay) 
■qn (strong decay) 
D,7r (strong decay) 
J/^Tvn (strong decay) 


A(1405) 
KKN 
DN 


1405 2 1/2- 
1920 4 1/2 1/2+ 
2790 2 1/2- 


qqs (L = 1) 


qqqsq 
qqqss (L = 1) 
qqqqc 


KN 
KKN 
DN 


20.5(R)-174(B) 
42(R) 
6.48(R) 


7rE (strong decay) 
Kir's, nr/N (strong decay) 
K + n~ir~ + p 


KNN 
DNN 


2352 2 1/2 0- 
3734 2 1/2 0- 


qqqqqs (L = 1) 


qqqqqq sq 
qqqqqq qc 


KNN 
DNN 


20.5(T)-174(T) 
6.48(T) 


AN (strong decay) 
K + tv~ + d, K + tt~tt~ + p + p 



TABLE II: Yields in one unit of central rapidity with oscillator frequencies uj — 550 MeV, uj a = 519 MeV, and uj c — 385 MeV. 





RHIC 


LHC 




2q/3q/6q 


4q/5q/8q 


Mol. 


Stat. 


2q/3q/6q 


4q/5q/8q 


Mol. 


Stat. 


/o(980) 


3.8, 0.73(ss) 


0.10 


13 


5.6 


10, 2.0 (ss) 


0.28 


36 


15 


a (980) 


11 


0.31 


40 


17 


31 


0.83 


1.1 x 10 2 


46 


U s (2317) 


1.3 x 10 -2 


2.1 x 10~ 3 


1.6 x 10" 2 


5.6 x 10~ 2 


8.7 x 10 -2 


1.4 x 10~ 2 


0.10 


0.35 


X(3872) 




4.0 x 10 -5 


7.8 x 10" 4 


2.9 x 10 -4 




6.6 x 10 -4 


1.3 x 10 -2 


4.7 x 10 -3 


A(1405) 


0.81 


0.11 


1.8-8.3 


1.7 


2.2 


0.29 


4.7-21 


4.2 


KKN 




0.019 


1.7 


0.28 




5.2 x 10~ 2 


4.2 


0.67 


DN 




2.9 x 10" 3 


4.6 x 10" 2 


1.0 x 10" 2 




2.0 x 10" 2 


0.28 


6.1 x 10" 2 


KNN 


5.0 x 10 -3 


5.1 x 10 -4 


0.011-0.24 


1.6 x 10 -2 


1.3 x 10 -2 


1.4 x 10 -3 


0.026 - 0.54 


3.7 x 10" 2 


DNN 




2.9 x 10" 5 


1.8 x 10" 3 


7.9 x 10" 5 




2.0 x 10" 4 


9.8 x 10" 3 


4.2 x 10" 4 



are used, and the temperature, chemical potentials are 
varied in a reasonable range. Moreover, the correlated 
uncertainties in the number of charm quark and its fu- 
gacity largely cancel out in the studied ratios. 

Our results also indicate that the yields of many multi- 
quark hadrons are large enough to be measurable in 
experiments. In particular, the heavy exotic hadrons 
containing charm or strange quarks can be produced at 
RHIC with appreciable abundance and even more so at 
LHC. Moreover, since the newly proposed states with 
charm quark are below the strong decay threshold, the 
background of their weak hadronic decays could be sub- 
stantially reduced through vertex reconstruction. Since 
the expected number of D° observed through the ver- 
tex detector is of order 10 5 per month at LHC, even the 
DNN states is definitely measurable. Therefore, rela- 
tivistic heavy ion collisions provide a good opportunity 
to search for multi-quark hadrons, and it may very well 
lead to the first observation of new multi-quark hadrons. 

In Fig. [I] we show the ratio Rh of the yields at RHIC 
calculated in the coalescence model 7V™ al to those of the 
statistical model Nf, tllt for the hadrons given in Table I. 



The grey zone within the range of 0.2 < Rh < 2 denotes 
the range of the ratios for normal hadrons with 2q and 
3g which are denoted by triangles inside the grey band. 
The ratios for the crypto-exotic hadrons with usual 2q/3q 
configurations also fall inside the grey band. The circles 
indicate the ratios obtained by assuming hadronic molec- 
ular configurations and are found to lie mostly above 
the normal band (Rh > 2). Moreover, we find that 
loosely bound extended molecules with larger size would 
be formed more abundantly. One typical example is 
A(1405). Using the previous relation between the bind- 
ing energy and the oscillator frequency oj, we find a small 
size for A(1405) (u = 174 MeV) and a ratio R h = 1.1. A 
coupled channel analysis 2§-3l[ gives, however, a larger 



(r ), leading thus to a larger Rh = 4.9. The patterns 
shown in Fig. [1] also holds for LHC as the freezeout con- 
ditions are similar to those at RHIC. 

As shown by diamonds in Fig. [T] the ratio Rh is below 
the normal band (Rh < 0.2) when a hadron has a com- 
pact multi-quark configuration. In particular, for light 
quark configurations, these ratios are order of magnitude 
smaller than those of normal hadrons or molecular con- 
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FIG. 1: (Color online) Ratio of hadron yields at RHIC in the 
coalescence model to those in the statistical model. 



figurations. This is consistent with the naive expectation 
that the probability to combine n-quarks into a compact 
region is suppressed as n increases. The tetraquark states 
of /o(980) and ao(980) are typical examples. This sup- 
pression also applies to 5q states in multi-quark hadrons 
(A(1405) and KKN) and the 8q state in KNN. On 
the other hand, the yield of hadrons at higher transverse 
momenta is expected to be enhanced if they have multi- 
quark configurations since quark coalescence enhances 
the ba ryon / m eson ratio at intermediate transverse mo- 
menta 15- 17 as observed in experiments [l|3, H3] ■ 

We conclude from the above discussions that the yield 
of a hadron in relativistic heavy ion collisions reflects 
its structure and thus can be used as a new method to 
discriminate the different pictures for the structures of 
multi-quark hadrons. As a specific example, we con- 
sider /o(980). So far STAR has a preliminary mea- 
surement of ,fo (980) /it and p / 71 " from which we find 
/o(980)/p° ~ 0.2 [32|. Using the statistical model predic- 
tion for the yield of p° = 42 leads to / (980) ~ 8. Com- 
paring this number to the numbers predicted for fo (980) 
in Table II, we find the data consistent with the KK pic- 
ture. Therefore, despite the quoted experimental error of 
around 50%, the STAR data can be taken as an evidence 
that the / (980) has a substantial KK component, and 
a pure tetraquark configuration can be ruled out for its 
structure. Such conclusion could not be reached from an- 
alyzing the data for / (980) — > 27 @,[33|. Because of the 
large error bars in the STAR data, further experimental 
effort is highly desirable for putting an end to this con- 
troversial issue. Similarly, efforts to measure the yields 
of other hadrons and newly proposed exotic states listed 
in Table I will provide new insights to a long standing 



challenge in hadronic physics. 
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